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ABSTRACT 

We have used high-speed spectroscopy of the rapidly oscillating Ap (roAp) star 
HD 99563 to study the pulsation amplitude and phase behaviour of elements in its 
stratified atmosphere over one 2.91-d rotation cycle. We identify spectral features 
related to patches in the surface distribution of chemical elements and study the 
pulsation amplitudes and phases as the patches move across the stellar disk. The vari- 
ations are consistent with a distorted nonradial dipole pulsation mode. We measure a 
1.6 kms~^ rotational variation in the mean radial velocities of Ha and argue that this 
is the first observation of Ha abundance spots caused by He settling through suppres- 
sion of convection by the magnetic field on an oblique rotator, in support of a prime 
theory for the excitation mechanism of roAp star pulsation. We demonstrate that 
HD 99563 is the second roAp star to show aspect dependence of blue-to-red running 
wave line profile variations in Nd ill spots. 
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1 INTRODUCTION 
1.1 The Ap stars 

Ap stars are chemically peculiar stars that range from early- 
B to early-F spectral types with the majority having de- 
tectable magnetic fields with strengths of a few 10^ to a 
few 10"^ G (iBvchkov. Bvchkova. Madeill2003l . according to 
whom 55 per cent have mean longitudinal fields greater 
than 400 G). The study of Ap stars with strong magnetic 
fields has ramifications for several other branches of as- 
trophysics. The interaction among strong magnetic fields, 
atomic diffusion and energy transfer in (or above) the upper 
atmospheres of non-degenerate stars has direct implications 
for observed stellar abundances and for the instability of 
stellar pulsations in /3 Cephei and sdB stars. The very differ- 
ent magnetic field strengths, orientations and geometries, as 
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well as the great variety of surface abundance distributions 
for Ap stars, therefore makes them particularly informative. 
The magnetic field geometry, determined from maps of lon- 
gitudinal magnetic field strength versus rotation phase, gen- 
erally agrees with an axisymmetric dipole, although strong 
deviations from this do o ccur, e.^. the quadrupole fiel d 
configuration of HD 37776 (|Thompson h Landstree^ll985l ). 
Typically, the magnetic axis is significantly inclined to the 
rotation axis, hence these stars are called oblique rotators. 

Many Ap stars are known to have surface spots, indi- 
rectly observable from photometric light variations or spec- 
troscopic line strength variations over a rotation cycle. These 
stars a re known as o ? CVn stars after the prototype of that 
name (jPyperl [19691 ). Their light and line strength varia- 
tions are a consequence of abundance spots that are of- 
ten associated with the magnetic poles. The rotation pe- 
riod of Ap stars is typically relatively long - of the order 
of several days to years, and even in some cases decades, 
as determined from their magnetic field variations. Slow 
rotation is a requirement for the atmospheric atomic dif- 
fusion that gives rise to the observed spectral peculiari- 
ties, since rapid rotation generates turbulent meridional cir- 
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culation that inhibits diffusion. Interestingly, atomic diffu- 
sion processes can give rise to observably stratified atmo- 
spheres in some Ap stars, and probably all rapidly oscillat- 
ing Ap (roAp) stars. Examples of studies of atmospheric 
abundance stratification in th e cool roAp star subgrou p 
are given bvlWade et all (l2QQlh: iRvabchikova et al.1 (|2002h : 
iRvabchikova. Leone Kochukhovl (|2005l ) . 

These studies are in agreement with each other. In 
general, Fe is concentrated by gravitational settling in 
the deeper observable layers in the photosphere around 
logrsooo ^ —0.5; the Pr and Nd line-forming layers are con- 
centrated by radiative levitation high in the atmosphere at 
logTsooo ^ —4, a level that is chromospheric in the Sun, 
and in the roAp stars is above the line-forming layer of the 
narrow core of the Ha line at —4 ^ logrsooo ^ —2, which is 
itself above the Fe line-forming layer. Many other elements 
show similar behaviour to these examples. 

Doppler imaging studies of Ap stars show that the spots 
associated with the magnetic poles are typically caused by 
abundance concentrations of rare earth elements. Overabun- 
dances globally are several order s of m agnitude compared 
to normal abundance stars (e.g., lAdelm an 1973), and may 
reach a factor of 10^ in the spots, as shown by the Doppler 
imaging study of the roAp star HR3831 by iKochukhovl 
(l2QQ4l ). and as shown in this paper for HD 99563. As a result 
of these strong abundance anomalies, line-blocking changes 
the temperature gradient at the magnetic poles of Ap stars; 
flux is redistributed usually from the blue to the red, so 
that rotational light variations in visible light occur in an- 
tiphase for short and long wavelength observations, along 
with obvious spectral line strength variations. Importantly, 
the rotational variations - whether magnetic, spectroscopic 
or photometric - allow the precise determination of the stel- 
lar rotation period. 

With the surface of some Ap stars resolvable by Doppler 
imaging, and atmospheric depth being resolved because of 
stratification, Ap stars - and particularly the roAp stars - 
offer the only opportunity for 3-D studies of stellar atmo- 
spheres for any stars other than the Sun. 



1.2 The roAp stars 

A subgroup of Ap stars, the rapidly oscillating Ap stars 
(roAp), exhibit rapid oscillations with periods of 5 — 21 min. 
Their frequencies are considerably higher than those typ- 
ical for 5 Set stars with whom they overlap in location 
inside the classical instability strip of the HR diagram. 
At present, only 40 roAp stars are known (see, e.g., 
iKurtz et al1l2006bl : IConzalez et al"]l2008l ). although several 
photometric sur veys have sear ched for rapid pulsation in 
Ap sta r s, such as Nelson Kre idl ( 1993 ); [Martinez Kurtzl 



Handler k Pa unzen (1999); Ashoka et ajj (|2000l ): 
IWeiss et al.1 (l2000l ): iDorokhova Dorokhov (20ol). For a 
(non-exhaustive) list of spectroscopic studies of roAp stars, 
see Kurtz, Elkin & Mathys (2006a). 

The roAp class characteristics are therefore poorly con- 
strained. As they share the same region of the HR dia- 
gram with the (photometrically established) non-oscillating 
Ap (noAp) stars, the selection of new roAp candidates is 
difficult and mostly le ads to null-results for bot h spectro- 
scopic studies (see, e.g. JPrevhammer et al ]l2008al : Freyham- 
mer, Elkin & Kurtz 2008b: lElkin et al.ll2008af r and for the 



photometric surveys listed above. The pulsations of roAp 
stars are described well by the oblique pulsator model, 
where the pulsation axis is assumed to be aligned with an 
oblique magnetic field axis (J<;urtz 1982, Saio 2005), or de- 
rived to be offset from both t he rotational and magnetic axes 
(iBigot k Dziembowskjll2002l ). 

iBalmforth et al.l (|200lh used this model and assumed 
that convection was suppressed locally in the magnetic po- 
lar regions. They showed that pulsation aligned with the 
magnetic poles is unstable to the axisymmetric, nonradial 
high-order modes observed in roAp stars, mainly through 
excitation by the opacity [n) mechanism acting in the hy- 
drogen ionization zones in these regions. They also consid- 
ered the effect of local He settUng, and concluded it is not 
the primary process responsible for the high-frequency pul- 
sation. Nevertheless, their model finds He to be deficient in 
the observable atmosphere at the magnetic poles. While He 
is not directly observable at visible wavelengths at the effec- 
tive temperatures of the roAp stars (6600 ^ Tgff ^ 8200 K) 
its gravitational settling implies that the magnetic poles are 
hydrogen- rich. We return to this point in Section [3.21 where 
we show evidence to support it. 

Statistical studies of the roAp stars can be intrigu- 
ing and informativ e, e.g. in the comparison of the roAp 
and noAp stars bv iHubrig et all (l2000l ). But with only 40 
such stars known the class characteristics are still not well- 
defined. Therefore, the greatest progress at present is to be 
made by in-depth stu dies of individual roAp stars, such as 
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that for HR3 831 fe.g. iKurtz et a 
and HR1217 (|Kurtz et al.ll2005al : ' 

This is reminiscent of the related problem in the pulsa- 
tional studies of 5 Set stars. Despite being the most abun- 
dant type of variable star, the major breakthroughs have 
been ob tained through ind ividual studies such as that of 
FGVir (iBreger et al.l l2005l ). For the ^Sct stars, as for the 
roAp stars, the question of when pulsation modes are ex- 
ited or not remains unanswered for physically similar stars 
- why do only some of these stars pulsate and what is the 
underlying physics? The many null results among roAp can- 
didate stars strongly suggest that either the pulsation prop- 
erties are not fully known (i.e. amplitude ranges or pulsation 
geometry) or the selection of candidates is too coarse. Indi- 
vidual studies of known roAp stars are therefore important 
(presently the only option at hand) to study the roAp prop- 
erties in detail. Some examples of enlightening studies in- 
clude that of 7 Equ that clearly demonstrated the trends for 
pulsational a mplitude with vertical stratification of chemi- 
cal elements (ISavanov, Malanushenko, Rvabchikoval 19991 : 
iKochukhov Rvabchikoval bOOll V Doppler imaging of 
HR3831 (Kochukhov 2 004, 2006), and the study of line 
bisector variability by iKurtz et al.l (l2006c) and a new 
type of upper at mospheric pulsation (|Kurtz et al.l l2006al : 
iKurtz et al]l2007l ). 



1.3 The roAp pulsations 

To first order, pulsations in many roAp stars are oblique ax- 
isymmetric dipole (/ = 1, m = 0) p modes. Nevertheless, it 
has been clear for decades that normal modes described by 
single spherical harmonics are insufficient to explain fully the 
mode geometry in these stars. Uniquely, the oblique modes 
of the roAp stars allow them to be viewed from varying as- 
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pect, giving detailed information about the mo de geometry. 
An ea rly example was the study of HR 1217 bv lKurtz et al.l 
(|l989l ) who showed the presence of alternating even and odd 
degree p modes, with the odd modes being consistent with 
oblique dipole modes, and the even modes not being consis- 
tent with either radial or quadrupole modes. This was con- 
firmed by the more extens i ve Wh ole Earth Telescope study 
of this star bv lKurtz et al.l (|2Q05al V Even for roAp stars with 
single pulsation modes, in the best observed cases - such as 
HR3831 - it can be shown that the single mode is a dis- 



torted oblique dipole mode (|Kurtz et alj|l997l : iKochukhovl 
I2OO6I ). 

Any deviation from spherical symmetry for a pulsat- 
ing star will give rise to distorted modes. For the roAp 
stars the dominant effect is the magnetic field, with contri- 
butions from the non-uniform abundance distributions and 
stellar rotation. The non-uniform abundance distributions 
are particularly important for spectroscopic st udies - such 
as the one in this paper on HD 99563, and that of lKochukhovl 
(|2006l l for HR3831 - where concentrations of ions in spots 
mean that the pulsation mode is non-uniformly sampled 
with the rotation. However, photometric studies that more 
uniformly sample the mode geometry with rotation also 
show the distortions of the modes. 

The effects of the magnetic field on the oscillations 
of roAp stars have been extensively studied theoret- 
ically, and there are ongoing investigations into this 
comp lex problem (jPziembowski Goode 19961: iBigot et al.[ 
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effect on the oscillations is expected to be small, the theo- 
retical result s do lead to expe c tation s of distorted modes. 
For example, ISaio Gautschvl (|2004h studied the infiuence 
of the magnetic field on the pulsation geometry and found 
that not only can a dipole magnetic field explain the 
suppression of the observationally absent 6 Set modes, but 
it also leads to insufficiency of a single spherical harmonic 
to describe the pulsations angular dependency; a series 
of spherical harmonics of different degrees / are required. 
Inclusion of magnetic fields in the study and interpretation 
of roAp pulsations is therefore necessar y and constitutes 
^a fo rmidable mathematical problem' (|Saio Gautschvl 
I2OO4 



[Sli3 (|2005l ) continued the study of the effects of a dipole 
magnetic field on roAp pulsations and showed again that the 
field stabilised low-order p modes (such as those excited in 
6 Set stars) for fields stronger than 1 kG, but that nonra- 
dial, high order distorted dipole, or quadrupole, p modes 
remain overstable and are most likely to be excited in roAp 
stars. This is in excellent agreement with the current obser- 
vations, as well as with the mounting observational and the- 
oretical evidence for observable pulsation nodes in and above 
the photosphere in the magneto- acoustic layer. Further, and 
most interestingly, Saio modelled the latitudinal amplitude 
dependence for the roAp case of HR383 1 and found reason - 
able agreement with the observations bv lKochukh^(|2004l ). 

Saio's non-adiabatic analysis predicts (his figure 10) a 
latitudinal dependence of amplitude that does not differ 
from a purely dipole mode photometrically, whereas spec- 
troscopy's ability to distinguish potentially between hori- 
zontal and radial pulsation components may detect the rel- 



atively smaller horizontal component when studying stel- 
lar spots close to (but not necessarily at) the magnetic 
poles. Saio's figure 10 demonstrates a reasonable agreement 
with observed latitudinal velocity amplitude variation for 
HR3831, but predicts higher amplitude toward the mag- 
netic axis than is observed. However, the formation layer 
of the Nd lines used by Kochukhov in his observat i ons is 
far above the boundary of Saio's models. iKochukhovl (|2006h 
also did not take into account the inhomogeneous Nd surface 
distribution. 

The oblique pulsator model describes how nonradial 
pulsation modes seen from varying aspect with the stel- 
lar rotation show amplitude and pha se modulation (see the 
overview bv iKurtz Martinez! l2000h that results in a fre- 
quency multiplet in an amplitude spectrum, centred on the 
pulsation frequency u and separated by exactly the rotation 
frequency z^rot • This exact splitt ing by the rotation freque ncy 
allowed Kurtz (1982; see also iKurtz Martinez! I2OOOI ) to 
rule out the possibility that the observed multiplets can be 
explained by rotationally perturbed m-modes (m being the 
azimuthal order of a mode), which is the typical interpreta- 
tion of frequency splitting observed in other (non-oblique) 
nonradially oscillating stars, as such rotational multiplets 
are split by (1 — Cni)^, where Q is the rotation frequency, 
Cni is a small constant dependent on the stellar structure, 
n is the overtone of the mode and / is its degree. 

For oblique pulsation the amplitude ratio of the mul- 
tiplet components contains information about the pulsa- 
tion geometry (assumed to be aligned with the magnetic 
geometry). In the case of a simple oblique dipole mode 
(A+i + A-i)/Ao = tan i tan /3, where i is the inclina- 
tion angle of the rotation axis, f3 is the obliquity of the 
pulsation/magnetic axis relative to the rotation axis, and 
Ao, A+i, A-i are the observed amplitudes of the frequency 
triplet. For modes that are not purely dipolar, or vary in 
structure with atmospheric depth in the stratified roAp stars 
atmospheres, this constraint can depend on the atmospheric 
depth of formati on of the spectral lines studied (see, e.g., 
lKochukhovll2006h . 

Despite the complex nature of roAp pulsations in the 
presence of magnetic fields and an oblique geometry, the 
observational data have rich potential for asteroseismic in- 
ference. Photometrically, multisite observations and space- 
borne observations allow the resolution of both the large and 
small frequency spacings and through them constrain the 
stellar dimensions and the interaction of the magnetic field 
and rotation with the pulsation modes. High-resolution time 
series spectroscopy allows the study of the pulsation geom- 
etry over the stellar surface (through line profile variability) 
and also vertically in the stellar atmosphere, through radial 
velocity measurements at different line bisectors, or through 
lines of elements loc ated at diffe rent depths in the strati- 
fied a tmosphere fe.g. jKurtz et al. 2006c; Rvab chikova et aP 
I2OO7I ). It is clear that the observations have begun to re- 
solve the magneto-acoustic nature of the pulsations in these 
stars high into the atmosphere to levels that are unobserv- 
able in any other star but the Sun. Radial nodes are re- 
solved in some stars such as HD 137949 (Mkrtichian, Hatzes 
& Kanaan 2003; Kurtz, Elkin Mathys 2 005c) and possibly 
HD 99563 (|Elkin. Kurtz Mathvsll2005l ) . Observed varia- 
tion of pulsation phase with atmospheric depth shows the 
presence of running waves, and there is even an indication 
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that s hock waves are observe d in the upper atmosphere of 
7Equ (|Shibahashi et al.ll2008l ), and by inference other roAp 
stars with similar hne profile variations. 

1.4 HD 99563 

HP 99563 was discov e red to be a roAp star by 
iDorokhova Dorokhovl (|l998h , pulsating in a 10.70-min 
mode. It is a visual binary with a 1.2 ma^ fainter secon dary 
at 1.79arcsec separation (|Fabricius Makarovl lioool ). A 
spectroscopic study by Elkin et al. (2005) confirmed the 
monoperiodic pulsation found from photometry, and found 
the pulsations to have extremely high radial velocity am- 
plitudes for a roAp star, such as 2.6 kms~^ for Ha and 
up to 5 kms~^ for Eu and Tm. For sake of clarity, we em- 
phasize that monoperiodic refers to a single (physical) pul- 
sation mode - e.g., the frequency quintuplet detected by 
[Handler et aD (|2006l ) (actually 5 components of a frequency 
septuplet) are all associated with only a single mode with 
pulsation frequency the rest of the quintuplet describes 
the amplitude and phase modulation of the mode. 

Elkin et al. (2005) measured a projected rotation ve- 
locity of vsini = 28.5 ± 1.1 kms~^ from a high-resolution 
spectrum, while from narrow-band Stromgren photometric 
indices they estimated Teff = 7700 K and log^ = 4.2 (cgs). 
Furthermore, they studied variations in amplitude and phase 
for line bisectors, and from element to element. The demon- 
strably high amplitudes (and thus high signal-to- noise, S/N) 
and rich and varied information content from using lines of 
different elements or with line depth, show HD 99563 to be 
a p articular interes t ing ro Ap star for detailed study. 

[Handler et all (|2006h carried out a photometric mul- 
tisite campaign and found an equally spaced frequency 
quintuplet in HD 99563, for the pulsation frequency u = 
1.5576539 ± 0.0000007 mHz. They determined the rotation 
period to be Prot = 2.91179 ± 0.00007 d and a time of 
pulsation amplitude maximum to be HJD 245 2031.29627. 
The photometry showed, through application of the a xisym- 
metric spherical harmonic decomposition method of iKurt j 
(|l992l \ that the mode of HD 99563 is dominated by a dipole 
component, with some contribution of an / = 3 component 
(i.e., the mode is a distorted dipole). 

Handler and collaborators combined photometry, a Hip- 
parcos trigonometric parallax and the temperature and sur- 
face gravity estimate by Elkin et al. (2005) to constrain the 
basic parameters of HD 99563. Under the assumption that 
the secondary star in the visual binary is a physical binary 
component and that the wide orbit indicates the stars have 
evolved separately, the estimated parameters of both stars 
must fit two stellar evolutionary models of the same age. By 
adopting Teff = 7900 di 300 K and log^ = 4.152, they found 
for the primary logL/L© = 2.03, R/Kq — 2.38, an inclina- 
tion angle of the stellar rotation axis i — 43? 6 zb 2?1, and a 
magnetic obliquity /3 = 86?4dz0?3. This makes HD 99563 es- 
pecially interesting because of its high pulsation amplitude, 
its rotation that is fast enough for Doppler imaging, and 
its (for a roAp star) rare favourable orientation where both 
ma gnetic poles become visible throughout a rotation cycle. 

iHubrig et all (l2004l ) were the first to measure the mag- 
netic field of HD 99563, finding a longitudinal field of — 688di 
145 G on HJD 245 2494.483, or rotation phase 0.08 with re- 
spect to Handler et al.'s pulsation ephemeris. iHubrig et al.l 



(|2006l ) published another two measurements at rotation 
phases 0.07 and 0.09 showing -235 ± 73 and 670 ± 37 G, 
respectively. Further, using the Russian 6-m telescope at the 
Special Astrophysical O bservatory of the R ussian Academy 
of Sciences fSAORAS). lElkin et al.1 (l2008bl ) obtained 6 ad- 
ditional magnetic longitudinal field measurements near ro- 
tation phases 0.25, 0.5 and 0.75. With all 9 available mea- 
surements, these authors used the well-known relation for a 
dipolar magnetic field: 

B\ oc Bp cos a (1) 
where 

cosa = cosi cos;^ + sini siuy^cos Qt, (2) 

a is the angle between the magnetic pole and the line-of- 
sight, i is the rotational inclination, (3 is the angle between 
the rotation axis and the magnetic axis, is the rotation 
frequency, B\ is the longitudinal magnetic field strength and 
Bp is the polar field strength. 

As can be seen from Eq.[2l the mean magnetic field 
strength {B\) oc cos i cos /3 and the amplitude of the mag- 
netic field variations Ab^ oc sin i sin /3. Thus tan i tan /3 = 
Abi/ {B\), from which P can be constrained when i is known. 
Using the measured valu es for HD 99563, (B i) = 21 di 58 G 
and Ab, = 701 ± 114 G, lElkin et al.1 (l2008bl ) find an incli- 
nation of i = 43? 5, and a magnetic obliquity P = 88° (84- 
90 ° for a la range), wh ich are similar to those determined 
bv lHandler et alT(|2006l V The magnetic minimum coincides 
with the moment of maximum modulated pulsation ampli- 
tude ((frot — hereafter), and similarly for the maximum 
half a rotation phase later {(frot — 0.5). This is consistent 
with both magnetic poles being viewed from a similar angle 
and coming into view for almost the same duration (requir- 
ing that either i or /3 be close to 90°), and with the pulsation 
axis and magnet i c axis being aligned. 

iKochukhovl (|2004l ) used Doppler imaging to construct 
the first tomographic maps of pulsation velocity, surface 
abundance and magnetic field map for the roAp star 
HR3831. With observations obtained over a wide range in 
rotation phase for HD 99563, it will be possible to do this 
with higher precision for this star because of its very high 
amplitudes, and it will be possible to look simultaneously at 
the additional dimension of the atmospheric depth depen- 
dence of the velocities. Kochukhov, as a first approximation, 
treated the analysis of HR3831 as if the atmosphere has a 
single pulsation layer. Observations of HD 99563 and other 
roAp stars show the situation to be much more complex 
than this. Disentangling the vertical and horizontal pulsa- 
tion characteristics is a challenge, but one with rich rewards 
for the fields of asteroseismology, pulsation theory, magne- 
tohydrodynamics, stellar atmosphere theory and atomic dif- 
fusion. 

HD 99563 is one of only a few known roAp cases where 
tomographic mapping through Doppler imaging is possible, 
based on rotation rate and the orientation of the magnetic 
field. Therefore, observations over three nights (most of the 
2.91-d rotation period) were scheduled from two sites si- 
multaneously - Subaru on Mauna Kea and VLT on Cerro 
Paranal - to attempt to collect spectra at high resolution 
in both time and dispersion to obtain a first picture of this 
star's photosphere, and to determine precisely its physical 
parameters. Theoretically, the scheduled observations would 
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Table 1. Journal of VLT observations indicating 
nightly observing date (UT), rotation phase with ref- 
erence to a time of spectroscopic pulsation maximum, 
HJD 245 4171.43328, Heliocentric Julian Date (HJD) 
range, typical exposure time, number of collected spectra, 
and mean S/N. The S/N ratios are based on the random 
noise in the continuum of consecutive spectra, or for sin- 
gle spectra based on predicted ratios. The 2004 data were 
re-reduced in same manner as the new data, and included 
in the analysis as they supplement the covered rotation 
phases. An additional 35 spectra were obtained with the 
HDS spectrograph at the Subaru telescope on 2007 March 
12 and March 13 and combined into two average spectra 
at rotation phases cprot = 0.23 and 0.56. 



Date (UT) Rotation HJD range texp n S/N 

phase cprot (—245 0000) (s) Asioo A6400 

2004 Mar 06 0.937-0.966 3070.60-3070.69 40 111 120 95 

2007 Mar 12 0.021-0.150 4171.53-4171.91 40 518 75 55 

2007 Mar 13 0.365-0.494 4172.52-4172.90 40/80 417 85 60 

2007 Mar 14 0.708-0.837 4173.52-4173.90 65 372 110 75 



have given us 70 per cent coverage of the rotation cycle, 
hence of the visible stellar surface. But due to bad weather, 
only one site contributed significantly to the pulsation study, 
and only some 40 per cent rotational coverage was attained. 

In this first of a series of papers on HD 99563, we are 
studying the horizontal pulsation geometry with the stellar 
rotation. A Doppler imaging study of the surface abundance 
distribution of many ions and detailed modelling of the hor- 
izontal and vertical pulsation geometry is in preparation. 



2 OBSERVATIONS AND DATA REDUCTION 

Time-series spectroscopy observations were collected on 
Cerro Paranal with the VLT UV- Visual Echelle Spectro- 
graph (UVES) during three nights in 2007 March; a jour- 
nal of observations is given in Table[T] UVES is an echelle 
spectrograph and we used the settings for the wavelength 
range 4970 - 7010 A (with a 60-A wide gap near 6000 A) at 
R = 110 000. The UVES instrument was used with an im- 
age sheer and a 0.3arcsec slit. Typical exposure times were 
40 and 65s, and 80s at high airmass on the second night. 
With ^ 26 s readout and overhead time, a time resolution 
of 66 — 106 s was obtained. 

Table[T] lists the spectra collected each night. With the 
2.91-d rotation period of HD 99563, the 3 night run was 
scheduled to cover most of a rotation cycle of HD 99563 
(rotational phases 0.02 — 0.84) and was coordinated with 
a run using the High Dispersion Spectrograph (HDS) on the 
8.2-m Subaru telescope at Mauna Kea to obtain the rota- 
tional phases coinciding with daytime on Paranal. However, 
humid weather conditions at Mauna Kea resulted in only 
35 sporadic observations that are not included in the fre- 
quency analysis in Sect. 13. ll but were instead combined into 
two average spectra at rotation phases 0.23 and 0.56 and 
used in Doppler imaging to produce surface maps for chem- 
ical elements, only two of which are shown in this paper 
(see Fig. [8] and Fig. 1131 below): a full Doppler imaging anal- 
ysis will be presented in a future publication. For details 



on the Doppler imagi ng programme and production of the 
tomographic maps see ISavanov Strassmeierl (|2005h . 

The spectra were reduced to 1-D with the pipelines pro- 
vided by ESO with the standard calibration data (bias, flat- 
field and Thorium- Argon wavelength reference spectra) . The 
S/N was 27 per cent lower than predicted with ESO's on- 
line exposure time estimation tool. The observing conditions 
were in general good, all 3 nights were clear and dry, with 
typical seeing of 0.5 — 0.8arcsec, except for the first night 
where seeing was 0.7 — 1.1 arcsec during the first half of the 
night and then deteriorated to 0.8 — 2.4 arcsec thereafter. 

Special effort was therefore made to optimise the reduc- 
tion. We used the UVES pipeline 2.9.7 (based on MIDAS) 
and applied minor additional corrections: We smoothed out 
a CCD bias pattern for every 3.5 columns (RED upper 
CCD), inserted 2 columns missed during CCD readout (with 
minor effects on the 2-D wavelength calibration), and fixed 
the rebinning step sizes of the 1-D spectra to 0.01787 A 
(RED upper CCD) and 0.01526 A (RED lower CCD). To 
eliminate wavelength drifts, wavelength calibrations were 
made to the nearest Th-Ar reference spectrum (obtained 
every 3 hours). 

Continuum rectification of the 1-D spectra was per- 
formed iteratively: individual spectra were normalised to 
a run master from the most stable night (starting on 
2007 March 12) using a 4th-order polynomial, then these 
were normalised to master spectra of individual nights. Af- 
ter inspection, subsets of the resulting spectra were cor- 
rected for remaining slopes. The nightly masters were rec- 
tified through iterative normalisation using spline fitting to 
pseudo- continuum windows identified from comparison to a 
model spectrum. Features due to cosmic rays and blemishes 
were removed with a semi-automatic cutting procedure. 
Barycentric corrections, which changed by about 1.0 kms~^ 
over a single night, were registered and applied in the pe- 
riod analysis, but not in the wavelength calibration and final 
rebinned 1-D spectra. 

Tab le[T] gives the estimated signal-to-noise ratio (S/N) 
of the spectra based on random noise measured in differ- 
ence spectra of pair-wise subsequent spectra. Because the 
2007 data had 27 per cent lower S/N than predicted, and 
also compared to the 2004 observations for the same expo- 
sure time, exposure times were increased to 60 s most of the 
two last nights in 2007. This difference is not due to weather, 
as large seeing variations are found to have only small in- 
fluence on the S/N (as a consequence of the image sheer). 
Comparing the 2-D spectra of the two years confirms that 
fewer photons were collected in 2007 than in 2004. About 
10 per cent variation in S/N is due to the differences in air- 
mass, since in 2007 we observed HD 99563 during full nights. 
Part of the reason for the lower S/N is from centring of the 
sheer on the object, but alone this cannot explain the loss. 
The spectra from the upper CCD, redward of 6000 A, show 
a systematic noise pattern of undulations with ^3.4 cycles 
per A, for stretches of ~30 A occurring every ~67 A, starting 
at 6070 A. This coincides with the first half of each Echelle 
order. As HD 99563 rotates relatively quickly, resulting in 
broad lines, we smoothed 1-D spectra of the upper CCD 
over 3 — 4 pixels, which strongly suppressed the artefact 
pattern. Comparing our frequency analyses for lines from 
the lower and upper spectra gave no significant differences. 
We suspect the origin of this feature is a sort of amplifica- 
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tion of the aforementioned 3.5 column pattern at the order 
merging regions. During the last night a 1 h gap in the time 
series occurred when the image acquisition system failed to 
save the observations. 



3 DATA ANALYSIS AND RESULTS 
3.1 Line analysis 

For the purpose of spectral line identification, a syn- 
thetic comparis on spectrum was produced with SYNTH 
(jPiskunovUlOoi ) using a Kurucz stellar atmosphere model. 
Atomic line data were taken from th e Vienna Atomic Line 
Database fVALD. lKupka et allfToOOl ) for ions with increased 
abundances, mainly for Nd, Pr, Sr, Cr and Eu. Other sources 
used for line data were the atomic database NIST (National 
Institute of Standards and Technology )[3 and the Database 
on Rare Earth Elements at Mo ns University (dreaMI; 
iBiemont. Palmer i h QuinetlllQQOl ) through its implementa- 
tion in the VALD. The line selection and line analysis is fully 
based on lines identified in table 1 of Elkin et al. (2005). 

3.1.1 Radial velocity shifts 

Because of the non-uniform surface distributions and strati- 
fication of ions in the atmospheres of Ap stars, high-overtone 
mode pulsation amplitude and phase vary from element to 
element, and even within individual lines as a function of 
line depth (atmospheric height). Line profiles are variable 
because of nonradial pulsation, non-uniform abundance of 
elements over the surface because of spots, and non-uniform 
abundance with depth because of stratification. It is our in- 
tention with this and future studies of the data presented 
here to disentangle these contributions to extract 3-D radial 
velocity maps of the stellar pulsation geometry. 

IDL tools for line measurements and analyses were used 
as described in lFrevhammer et al.l (|2008al V Precise radial ve- 
locity shifts were measured for many spectral lines with the 
centre-of-gravity (CoG) method and by fitting with Gaus- 
sian profiles. We made our measurements with respect to the 
continuum - normalised to unity - in the selected regions of 
the measured lines, which gave smaller amplitudes for the 
CoG procedure. This method, however, was preferred for 
pulsation studies because of its better ability to deal with 
blended line profiles. For strong and isolated lines, the two 
methods were comparable. However, in some cases Gaus- 
sian fitting was used in automated measurements as it is 
more stable when the whole line profile shifts considerably 
due to, e.g., rotational mean radial velocity modulation (see 
Sect. 3.2) when one of the line wings goes outside the pre- 
defined wavelength window. The main difference of the two 
methods is their estimates of pulsation amplitudes. The CoG 
method is very sensitive to continuum placement and more 
sensitive to subtle, intrinsic line profile variations of rela- 
tively weak, asymmetric lines. 

The radial velocity measurements were automated with 
a wavelength window defined for each line in a high S/N 
average spectrum. That window was shifted in wavelength 

^ http://physics.nist.gov 

^ http:/ /w3. umh.ac.be/~astro/dream.shtml 



to account for the barycentric correction for each individ- 
ual spectrum, then the profile was fitted with both a Gaus- 
sian and by CoG fitting; finally, barycentric corrections were 
added to the resulting radial velocities. The two radial ve- 
locity series were then compared. 

Because of the occurrence of substructure in absorption 
features of elements with inhomogeneous surface distribu- 
tions, we used two approaches to study line variability: 1) 
to study an element's integrated radial velocity behaviour, 
we used the full line profile, whereas 2) for line profiles with 
sub-features caused by inhomogeneous (spotted) surface dis- 
tributions, we measured each substructure in a line and fol- 
lowed it over the half a rotation cycle for which it was visible. 
The line windows for the 'spot' features (two for Nd and 3 
for Eu, see below in the text and Figs[Tn]and[T4]for clarifica- 
tion) were defined in 1 hr sub-series of spectra (for a total of 
32 time bins), rather than one window for all 1418 spectra. 
The former approach was useful for fitting multiplet terms 
(Sect. 3.2), while the latter approach was particularly useful 
for studying pulsation amplitude and phase variations with 
rotation phase. 



3.1.2 Frequency analyses 

Frequency analyses were performed u sing a Discrete Fourier 
Transform programm e (iKurtd [TosBI ) and the PERIOD04 
(jLenz k Bregerll2005l ) programme. The former fits cosines 
while the latter uses sines, so all phases determined with 
PERIOD04 were recalculated to match the cosine fitting 
in units of radians. The noise, a, of fitted amplitudes and 
phases (see, e.g., Tables[2]and|4| was determined frora least- 
squares fitting to the data following iDeemin 3 (119751 ): it is 
the standard deviation of one measurement with respect to 
the fitted function. 

Trends in the radial velocities with the rotation were 
strong for all lines, even Ha and Fe (Fig.[T]), indicating that 
HD 99563 is very spotted. Linear and non- linear trends were 
fitted and removed with linear or polynomial least-squares 
fitting, and in a few cases also by prewhitening low frequency 
peaks from the radial velocity series. The origin of these 
trends are discussed further in Sect. 3.2, but most of the 
contributions are stellar; with the 0.3-arcsec slit, no image 
sheer and seeing conditions of 0.9 — 1. 4arcsec, the centrin g 
error for UVES is only 50 - 100 ms"^ (iBouchv et al.ll2004l ). 
Furthermore, a 1-mbar change in pressure may induce drifts 
of 90 ms~^. Therefore, instrumental shifts were less than 
^300 ms~^ per night, depending on seeing and pressure. 

The telluric lines in Fig.[T]show that the wavelength cal- 
ibration sometimes has small offsets (up to 600 ms""*^) be- 
tween reductions using different wavelength reference spec- 
tra. These are artefacts from not interpolating over all ref- 
erence spectra for a night, an option not available with the 
UVES pipeline, and then rebinning to same- width wave- 
length bins. However, as the time scale of these offsets is 
very different from that of the pulsations (3hr as opposed 
to llmin), we did not correct this in the analysis. Telluric 
fines only provide accurate and precise velocity references 
within a few hundred ms~^, as they are infiuenced by fast 
and changing wind speeds in the high layers of the Earth's 
atmosphere where they are formed. 
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Figure 1. Rotational radial velocity trends for different elements: 
telluric lines (2 lines combined), Mg (two Mgl lines combined), 
Fe (combined from 4 Fel and Fell lines), and the hydrogen core 
(prewhitened for pulsations). Barycentric corrections were applied 
to the radial velocities of stellar lines and are relative to those of 
the telluric lines. Radial velocity offsets were added to each series' 
ordinates for clarity. If the radial velocity variation with the rota- 
tion period was the result of reflex motion of HD 99563 in a binary 
orbit with an unseen close companion, the amplitudes and trends 
of relatively homogeneously distributed elements, such as iron and 
hydrogen, should be the same; they are not. The radial velocities 
for the iron measurements that deviate most from those for Ho; 
at phase 0.94 are from 2004, but there is also signiflcant deviation 
at phases 0.02 — 0.17 in the 2007 data. Mg shows signiflcantly dif- 
ferent trends from Fe and H, but this is most probably caused by 
a patchy surface distribution. Because of the non-uniform abun- 
dance distributions, orbital radial velocity shifts of HD 99563 can 
only be ruled out for amplitudes signiflcantly larger than those 
seen here. 



3.2 Hydrogen core 

iHandler et alJ (|2006l ) detected 10 frequencies in their time 
series multisite campaign B-hand photometry of HD 99563. 
These are comprised of (Table[2| the pulsation frequency 
u = 1.557653 ± 0.000007 mHz, the rotational split oblique 
dipole frequency quintuplet: zyiz^rot and z^zb3z^rot, and other 
harmonics or combination frequencies. These authors iden- 
tified the rotation frequency as Urot = 0.0039751 mHz and 
gave the photometric ephemeris (in days) with respect to the 
time of photometric modulated pulsation amplitude maxi- 
mum as: 

t(max) = HJD 245 2031.29627 + 2.91179 ± 0.00007 E, (3) 

where t(max) gives the time in HJD corresponding to the 
given number of epochs E elapsed since the reference point 
in time. The decimal part of E corresponds to the rotation 
phase in the case of HD 99563. The pulsation frequency 
ly itself has the second lowest amplitude of the quintuplet. 
See Handler et al.'s figure 5 for a schematic view of their 
frequency quintuplet. 

The Ha-core line-forming region extends over a con- 
siderable vertical region of the chemically stratified photo- 
sphere, in contrast to other elements. Radial velocity mea- 
surements of this line therefore compare well to the photom- 
etry, and as Ha for this star has a high pulsation amplitude 
in excess of 2 kms~^, we used this line for the initial fre- 



quency analysis. Because of the limited accuracy of the ro- 
tational period in Eq.[3l a reference epoch close to our 2007 
observations was chosen at the moment of maximum of the 
modulated pulsation amplitude of Ha: 

t(max) = HJD 245 4171.43328 + 2.91179 (±0.00007) E, (4) 

while the v a lue of the pulsation frequency u by 
[Handler et aD (l2006l ) was revised within its error to v — 
1.5576530 mHz in order to tie the new spectra together with 
the 2004 data set (see Sect. [1211]) • 

3.2.1 Trends of the radial velocity shifts 

The Ha profile is shown in Fig. [2] The insert shows the 
Une core, with the region measured indicated with a bar. 
Radial velocities, measured through Gaussian fitting to the 
line profile cores, are shown in Fig. (3] phased with the ro- 
tation period. These reveal a double wave sinusoidal vari- 
ability of the average radial velocities with an amplitude 
of 1.6 kms~^ and with a frequency twice that of the ro- 
tation frequency, 2zyrot = 0.0079 498 mHz. Using the fre- 
quency analysis tool PERIOD04 (jPenz Bregerll2005l ) to 
i mprove phases and am plitudes iteratively for the solution 
of [Handler et all (|2006l ) (for fixed frequencies, keeping only 
the significant ones), the 2z^rot modulation is found (Table[2]) 
to have an amplitude of 1600 ms~^. Single or double wave 
angular modulation of the mean radial velocity of stars is 
known from many other elements in spotted stars - partic- 
ularly the rare earth elements - but has not previously been 
detected for hydrogen. 

Comparison with telluric lines, constant within 
600 ms~^, excludes a non-stellar origin (Fig.[T]) for the Ha 
rotational radial velocity variations. Detailed comparison of 
averaged spectra in bins of 35 minutes show the slow pe- 
riodicity in Fig.[T] dominates in the blue wing and central 
part of the slightly asymmetric Ha core. A blend with a 
line of a rare earth element distributed mainly around the 
stellar poles cannot be completely excluded as the source of 
the radial velocity curve. However, we consider it unlikely, 
based on comparison to a synthetic spectrum, that such rel- 
atively weak lines (compared to the Ha core) can result in 
the amplitude observed. 

The variation in the mean radial velocity of the Ha line 
with rotation is shown in Fig.O with the pulsation variation 
present (in the top panel) and after prewhitening for the 
pulsational nonuplet (given in Table 2). There are several 
possible explanations for this rotational variation that we 
now examine. 



3.2.2 Variations in the mean radial velocity caused by line 
profile variations 

We modelled line profile variations for an oblique dipole pul- 
sation mode with a rotational inclination of i = 44° , pulsa- 
tion pole obliquity of /3 = 88°, an equatorial rotational ve- 
locity of 30kms~^ and a pulsational amplitude of 14kms~^ 
- values appropriate for HD 99563. The synthetic line pro- 
files were computed using an adapted version of t he line 
profile code implemented in FAMIAS (|Zimal l2006l . 2008). 
The underlying model approximates the velocity field by 
a spherical harmonic, assumes a Gaussian intrinsic profile 
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Figure 2. Average spectrum of the Ho; region. The insert shows 
the hne core, with the region measured for radial velocities 
marked with a bar. 
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Figure 3. Radial velocities for HD 99563 measured in the line 
core of Ha for the 2007 observations (black dots) and the 2004 
data (grey dots at phase 0.94). The modulated pulsation am- 
plitude is superposed on a rotational modulation of the average 
radial velocity (fitted with solid line). Rotation phases are relative 
to Eq.3] The top panel shows measured radial velocities, while the 
bottom panel shows rotational radial velocities that remain after 
removing the pulsational variations by prewhitening the pulsation 
frequencies (i.e., the pulsational frequency nonuplet) in Table[2] 
The broad band of data in the top panel is not caused by scatter; 
it is the pulsation compressed in time. Fig. [5] shows these data 
expanded in time where the pulsation cycles can be seen. 



that originates from a two-dimensional atmospheric layer, 
and is weighted with a fourth-order limb darkening law. We 
adapted the above mentioned code to include the effects of 
oblique pulsation on the line profile variations. 

We computed a first set of synthetic line profiles for 
a complete rotation cycle, assuming uniform element dis- 
tribution across the stellar surface. We then measured the 
radial velocities for the model data by computing the first 
moment of the line profile. The top panel of Fig.|4] shows 
that in this case the radial velocity amplitude is modulated 
whereas the rotational variation of the mean radial velocity 
of the model line is negligible. Therefore, the radial velocity 
variation with the rotation period for Ha is not a conse- 




0.5 
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Figure 4. Top panel: Radial velocities over the rotation period 
(black dots in electronic version) and its running mean (solid line 
at RV=0 kms"-*^) for model line profile variations of an oblique 
dipole pulsation mode and uniform element distribution with i = 
44°, f3 = 86°, vsini = SOkms"-*^ and a pulsation amplitude of 
14kms~-'^. The radial velocity amplitude is modulated due to the 
changing aspect angle of the pulsation mode, whereas the mean 
shows no significant variation. Abscissa is rotation phase. Bottom 
panel: Same as the top panel, but with a non-uniform element 
distribution that is higher near the magnetic poles, simulating 
H abundance spots. Both the radial velocity amplitude and its 
mean are modulated due to the stellar rotation in the same way 
as the observations of the Ha line. As for Fig.O the model data 
here are compressed in time so that the pulsations appear only 
as a broad band. If this figure were expanded in time (as in Fig. [5] 
for the Ho; data), then the pulsation cycles would be seen. The 
important features to be noted are the envelope of the data, which 
shows the amplitude modulation, and the mean radial velocity. 
Comparison of the lower panel with the Ho; data in the upper 
panel of Fig. [3] shows the plausibility of our hypothesis that the 
variation in mean radial velocity is caused by enhanced hydrogen 
spots at the magnetic poles. 



quence of line profile variations for a simple oblique dipolar 
pulsation mode. 



3.2.3 An unseen binary companion in a synchronous orbit 

An obvious possible explanation for the radial velocity mod- 
ulation seen in Fig. [3] is to hypothesise an unseen binary 
companion. In this case, the orbital period must be equal to 
the rotation period of HD 99563, and that is reasonable for 
a 2.91-d rotation period as synchronism is expected for such 
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a short orbital period. There is no secondary spectrum seen, 
thus the single-hned spectroscopic binary mass function is 

^^^^^^ = (v.sinif^ (5) 

(mi + 1712)^ 2tiG 

where the primary's projected orbital velocity visini — 
1600 ms~^ is from Table[2] The component masses mi and 
m2 are for the primary (pulsating) component and the less 
massive secondary component respectively, P is the orital 
period (assumed identical to the primary's rotation period). 
Assuming the orbital plane is the same as HD 99563 's equa- 
torial plane - a reasonable assumption for such a close sys- 
tem - gives the orbital inclination i to be the same a s 
the rotational inclination derived by iHandler et al.l (|2006l ) , 
i — 44°. With this inclination the very small mass function 
shows that mi >> m2. Taking the mas s of HD 99563 to be 
that estimated by iHandler et al.1 (|2006l ). M = 2M0, gives 
m2 = 0.025 M0, or about 26 Jupiter masses with a separa- 
tion between the components of about 11 R©. 

Thus one explanation of the radial velocity variation 
in Ha with rotation is a close hot Jupiter companion in a 
synchronous orbit. The current exoplanet listingj indicates 
that this is not unreasonable when compared to other known 
exoplanets. This, if true, is an exciting result, as no planet 
has been discovered orbiting a magnetic Ap star. 

If a planetary companion is the explanation for the Ha 
rotational radial velocity variation, then we expect all spec- 
tral lines to show the same radial velocity curve. Testing this 
for HD 99563 is complicated by the spotted nature of its at- 
mosphere, causing all spectral lines to vary with rotation, 
but with different amplitudes and phases depending on the 
abundance distributions. For example, other elements, such 
as Mg in Fig.[T] show different mean radial velocity curves 
(also when using telluric lines as reference for the wavelength 
calibration), whereas in the same figure, Fe shows some sim- 
ilarities with Ha. Because the rotational radial velocity am- 
plitudes of various elements caused by the spots are compa- 
rable to, or higher than, the rotational amplitude variation 
seen in Ha, a hot Jupiter remains a possible explanation, 
but not our preferred one. 

3.2.4 Spots of enhanced H at the magnetic poles 

We suggest this is the first detection of enriched H abun- 
dance spots near the magnetic poles, as is expected the- 
oretically from helium depletion by gravitational settling 
(jBalmforth et al.1 l200lh . Such spots give rise to extra ab- 
sorption superposed on the mean Ha core Une profile that 
is variable with the rotation of the star. When a hydrogen 
spot is in the centre of the disk, the radial velocity contri- 
bution from its extra absorption is zero. Then as it moves 
away from the centre the extra absorption is in the part of 
the disk that rotates away from us, giving rise to a red shift 
of the line. As the spot reaches the limb, the opposite spot is 
coming into view on the other side of the disk, increasingly 
cancelling out the redshift by a blueshift. At quadrature (or 
in this case, magnetic crossover), both spots cancel out the 
contribution to the mean profile and the radial velocity shift 
is zero. This is followed by a blue shift. This is what is seen 

^ http: / / exoplanets.org 



in Fig.[3l the phasing is correct; the double- wave is as ex- 
pected. 

To simulate the scenario of enriched element abundance 
spots near the magnetic poles, we followed the same proce- 
dure as in Sect. 3.2.2 above, but in this case with weighting 
of the local intrinsic line width W with respect to the lati- 
tude of the displacement field using the following arbitrary 
function: 

W{0) = Wo{l + cos^O) (6) 

The resulting shape of the radial velocity curve is displayed 
in the bottom panel of Fig.|4l The radial velocity amplitude 
as well as its mean value are modulated with an amplitude 
that is comparable to the observations of the Ha line. 

Thus our hypothesis that the observed Ha line core's 
radial velocity variation with rotation is caused by theoreti- 
cally expected hydrogen abundance spots is consistent with 
our simple model. Furthermore, and more importantly, we 
predict that if our hypothesis is correct, then hydrogen spots 
will be detectable with stellar rotation for other Ap stars 
with strong magnetic fields of similar configuration and ori- 
entation as HD 99563. We plan tests of this prediction. 

Of course, even if hydrogen was uniformly distributed, 
there would be some rotation variation in its equivalent 
width as a result of the cooler spots at the magnetic poles 
caused by supression of local convection by the vertical mag- 
netic fiux tubes and by the non-uniform distribution of other 
elements - primarily rare earth elements. However, in this 
case the contribution, if any at all, from the polar regions is 
a reduced absorption in the hydrogen lines (for stars cooler 
than A5-type stars) , which is quite opposite to our hydrogen 
spot explanation. This does not agree with what is seen in 
Fig-El In addition, the general shape and phase of the mod- 
ulated mean radial velocity curve from hydrogen spots must 
match that of other elements concentrated near the stellar 
poles. Nd is such an example (see the Doppler imaging map 
in Fig. [8]) and, except from different amplitude and that the 
lines become double at magnetic crossover, the phases of the 
modulation are similar (see Fig. llip . 

3.2.5 Detected frequencies 

The frequency resolution of our 2007 data alone, 4.8 / iHz, i s 
insufficient to improve the period bv [Handler et al.1 (|2006h . 
However, by revising the pulsation frequency in a test 
grid within the la error given by Handler et al., we could 
also include the 2004 data in the fit. This revision was done 
by recalculating the first seven frequencies by Handler et 
al. (Column 2 in Table[2| relatively to the values of the 
trial frequency zy, using fixed amplitudes. A small revision 
(+4.8 /iHz) of z/, chosen as the closest alias to their value 
of z/, gave a good fit to the frequency quintuplet and resid- 
uals of only 245 ms~^ per observation for all 2004 — 2007 
data and 10 ms~^ for the highest peaks in amplitude spec- 
trum of the residuals to the fit. The revised frequency is 
henceforth used as zy, but we emphasize that it is a conve- 
nient alias out of many choices that phases the 2004 data 
with the 2007 data, rather t han an improvement in preci- 
sion to the multisite result bv IHandler et al.l (|2006l V We note 
that combining the two data sets cannot improve the fre- 
quency determination more than with the 2007 data alone, 
due to the introduction of several ambiguous aliases that 
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Figure 5. Radial velocities for HD 99563 measured in the line core of Ha for the 2007 observations and the 2004 data (top). Each 
panel is for a separate night. The data were prewhitened for the rotational mean radial velocity modulations. Superposed in blue is the 
frequency solution in Table[2] using the distorted oblique dipole nonuplet splitting of u (1.5576530 mHz). 



are all good candidate frequencies. Yet, the 2004 data sup- 
plement the pulsation analysis at other rotation phases and 
were therefore usefully included. 

After prewhitening the full radial velocity series (all 
2004 — 2007 data) for the mean radial velocity modulation 
(at 2z/rot) and the first seven si gnificant frequencies in the 
solution bv [Handler et al.1 (l2006l ). we detected excess power 
in the residuals. By including nine terms spaced equally in 
frequency by the rotation frequency, we found all but one to 
be significant; we refer to these nine frequencies as a nonu- 
plet. Figs[3]and[5]show respectively the radial velocity series 
prewhitened for the nonuplet frequencies and for the rota- 
tional mean radial velocity modulation. The amplitude of 
hydrogen pulsational radial velocity variations ranges from 
a few hundred ms~^ to almost 4 kms~^. Fig. [6] illustrates 
the nonuplet frequencies, while Table[2] (columns 4 — 5) gives 
the corresponding amplitudes and phases determined from 
a linear least- squares fit. 




1.54 1.55 1.56 1.57 1.58 

Frequency (mHz) 

Figure 6. The frequency nonuplet detected in the Ha radial ve- 
locity series for HD 99563. Mean amplitude error is indicated with 
a bar (top, right). Note the multiplet structure is a result of the 
distorted dipole pulsator geometry, not rotationally perturbed m- 
modes. The frequency splitting is exactly the rotation frequency. 



3.2.6 Amplitude and phase variations 

Amplitude and phase variations were measured with radial 
velocities from Gaussian fits to the line core of Ha (Fig.[2|) 
in 32 bins of about 1 hr duration. The result is shown in 
Fig. [71 The superposed full line is for a rotationally modu- 
lated amplitude nonuplet pattern for z/, based on amplitudes 
and phases in Tabled The dashed line shows the same, but 
for a pure dipole triplet pattern only. T hese are based on the 
general (2/ + 1) multiplet expression in iKurtd ([l992) where 
a linear combination of axisymmetric spherical harmonics is 



used to describe the observed angular amplitude and phase 
variations. Amplitudes and phases of the multiplet terms 
were determined from a fit to the multiplet frequencies from 
full line profile radial velocity series for that particular line 
using all spectra. For symmetric lines such as Ha, we used 
Gaussian fitting. Fig. [7] demonstrates graphically what Ta- 
ble^ shows numerically: The pulsation mode is basically an 
obhque dipole mode, but significant distortion from a sim- 
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Table 2. Amplitudes and phases determined from a least- 
squares fit of the listed frequencies to the Ha radial ve- 
locity series. Column 3 gives the frequencies calculated 
relative to a value of u slightly revised from that of 
(jHandler et al.ll2006h (to tie together the 2004 and 2007 
data sets) and to i/rot- The frequency nonuplet was re- 
quired to have components separated by exactly the rota- 
tion frequency. A 2 cd~^low frequency periodicity that is 
believed to be an artefact was removed. Errors on ly and 
Urot from (|Handler et all l2006l ) are cr^ = 0.000007 mHz 
and cTi^rot = 0.0000001 mHz. Prewhitening for this solution 
leaves point-to-point scatter in the residuals of 212 ms~^. 
Gaussian fitting to the profile was used. The zero point 
for the time scale to = HJD245 4171.43328 was selected 
such that the phases u — Urot and u + Urot are equal; 
i.e. this is the time of pulsation amplitude maximum, as 
given in Eq.|31 
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u 


1.5576530 


1 


5576482 


156 ± 


31 


-2 


56 ±0 


20 


5.1 




1.5536779 


1 


5536733 


1726 ± 


32 


-2 


17 ±0 


02 


54.3 




1.5616281 


1 


5616231 


1303 ± 


32 


-2 


18 ±0 


02 


41.0 


i/-2l/rot 




1 


5496984 


121 ± 


22 


-2 


11 ±0 


18 


5.5 


Z/+2l/rot 




1 


5655980 


80 ± 


22 





42 ±0 


27 


3.7 




1.5457276 


1 


5457235 


449 ± 


20 


-2 


63 ±0 


04 


22.2 


Z/+3z/rot 


1.5695784 


1 


5695729 


105 ± 


20 


1 


25 ±0 


19 


5.3 






1 


5417486 


35 ± 


19 


1 


56 ±0 


55 


1.8 






1 


5735478 


123 ± 


19 


-1 


87 ±0 


16 


6.4 


2v 


3.1153060 


3 


1152963 


129 ± 


8 





95 ±0 


06 


15.9 


i^rot 


0.0039749 






53 ± 


11 


-3 


14 ±0 


26 


4.9 


2z/rot 


0.0079498 






1600 ± 


13 


-1 


90 ±0 


01 127.0 


2l/-2z/rot 


3.1073557 


3 


1073465 


72 ± 


8 





48 ±0 


11 


8.8 



pie dipole geometry is required to explain the radial velocity 
amplitude and phase modulation with rotation for hydrogen. 

3.3 Neodymium 

This element is very inhomogeneously distributed on the 
stellar surface. The Doppler imaging map in Fig. [8] shows 
two concentrated spots, co-located with the magnetic poles 
as defined at rotation phases 0.0 and 0.5 by Handler et al. 
In the analysis of Nd, we will first study the frequency con- 
tent of the Nd lines using the full line profiles, and then 
study rotational amplitude and phase variations for these 
two spots. For practical purposes (see, e.g.. Fig. 1101 below), 
we identify these spots as the A (rotational phases 0.2 — 0.8) 
and B (rotational phases 0.7 — 1.3) regions. 

S.S.I Observed line profile variations 

We find significant partial doubling of the Nd lines during 
the magnetic crossover phase, which allows us to study the 
pulsation of two discrete regions of the stellar surface, except 
for rotation phases near pulsation maximums when only one 
or the other spot is visible. Practically, a direct implication 
of this is that the pulsation amplitude does not go to zero 
for this element at magnetic quadrature, as it does for Ha. 
However, as the pulsations in these spots are in antiphase, 
they would cancel out if we could not resolve their individual 
line contributions. This concentration of certain ions, such 
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Figure 7. Amplitude and phase variations from Gaussian fits to 
the line core of Ho;. The superposed model in the full line is for 
a rotationally modulated amplitude nonuplet pattern for u (see 
text). For comparison, a pure dipole fit to the frequency triplet 
only is shown with the dashed (blue) line. The lack of spectra at 
^rot =0.25 gives the fit too much freedom here. 



as those of Nd, in spots on HD 99563 (and some other roAp 
stars) allows us to study the pulsation behaviour of only a 
part of the visible hemisphere of the star from varying as- 
pect, an opportunity provided by no other kind of pulsating 
star. Ultimately, this may allow for full 3-D mapping of the 
pulsation velocity field, with the disentangling of horizon- 
tal and vertical velocities - important for comparison with 
theoretical predictions (e.g., 'Saio"2005'). 

To evaluate the frequency content in the Nd lines, radial 
velocities were measured with Gaussian fits to the symmet- 
ric Nd lines Nd ill 5102.41 A and Nd ill 6327.24 A, using the 
full line widths, and they were then combined. Spectra near 
rotation phase 0.75 were excluded as that is where the A 
and B pulsation phases become too mixed, or where one of 
the components strongly dominates the other. 

The radial velocity series for each spot were each de- 
trended with a linear or polynomial fit prior to the fre- 
quency analysis (an example of this detrending is shown in 
Fig. [9}. That is justified as the average radial velocities of the 
two spots shift quasi- linearly with rotation phase, as seen in 
Fig.Hni A few low frequencies in the range 2 — 5 cd~^were 
also included in the least-squares fit to get better error es- 
timates. For the fixed nonuplet frequency solution, as for 
Ha given in TableO we find all components to be signifi- 
cant, including ly — Aurot- The residuals show a flat amplitude 
spectrum with no indication of other pulsation frequencies. 
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Figure 10. Line region in the vicinities of Nd ill 5102.41 A (left) 
and Nd III 6145.07 A (right) shown for all spectra sorted accord- 
ing to rotation phase of HD 99563. The spectra are normalised 
to unity in the continua (white colour) and absorption features 
appear dark or black. Four phase intervals not covered by our 
high speed spectroscopic observations appear with a grey-scale 
corresponding to 0.95. Two line features relate to the dominant 
surface patches of high Nd concentration. These are referred to 
as the A (rotational phases 0.2 — 0.8) and B (rotational phases 
0.7 - 1.3) regions. 



Figure 8. Doppler imaging map for Nd ill 6145.07 A. The abun- 
dances are on a scale of log A^h = 12; the sol ar abundance of Nd is 
1.41 on this scale tA splund, Grevesse, Sa uval 200^. Thus the 
highest overabundances of Ndlll in the spots rise to 10^-^ times 
solar at the centres of the two spots associated with the magnetic 
poles. These are evident at rotation phases 0.0 and 0.5, and they 
lie near the stellar rotational equator (latitude 0°) as a conse- 
quence of the magnetic obliquity being near 90° {(3 = 88°). The 
white colour below latitude —44° represents the part of the star 
that is never visible because of the i = 44° rotational inclination. 



E 6 



E 6 





1.6 

frequency (mHz) 



De-trended 
RV series 



[.. .J 




1.4 1.6 

frequency (mHz) 







2 3 

frequency (mHz) 



Figure 9. The effect of detrending the 'non-binned' radial ve- 
locity curve for Nd ill 5102.41 A. Top panel shows the amplitude 
spectrum of the raw radial velocity series (with zoom-in panel 
of the region with the nonuplet pulsational frequency content), 
while the bottom panel shows the amplitude spectrum for the 
detrended radial velocity series. 



Table 3. Amplitudes and phases fitted to combined radial 
velocities of two neodymium lines. The zero point from 
Ha in Eq.|3] was used and the different phases for the 
two main sidelobes ly^Urot correspond to different zero 
points for times of maximum amplitude modulation for 
H and Nd. In fact, this difference in time is 0.19464 d 
with the Nd maximum coming before the H maximum. 
This occurs 0.067 earlier in rotation phase, or 24 degrees 
earlier in rotational longitude - in good agreement with 
the Nd abundance map, Fig.El 

Id. Freq. Amplitude Phase S/N 

(mHz) (ms-i) (rad) 



Z/-2z/rot 
Z/+2l/rot 
Z/ — 3l/rot 
Z/+3l/rot 

2v 

2i/rot 
2z/-2z/rot 



1.5576482 
1.5536733 
1.5616231 
1.5496984 
1.5655980 
1.5457235 
1.5695729 
1.5417486 
1.5735478 
3.1152963 
0.0039749 
0.0079498 
3.1073465 



328 ± 64 
1131 ±71 
531 ± 72 
445 ± 50 
456 ± 50 
305 ± 43 
192 ± 43 
230 ± 43 
171 ± 43 
77± 18 
134 ± 31 
141 ±31 
158 ± 18 



-0.28 ± 
-2.87 ± 
2.57 ± 
1.09 ± 
-1.54 ± 
-1.19 ± 
-0.68 ± 
-3.04 ± 
2.02 ± 
0.68 ± 
-0.63 ± 
1.69 ± 
1.13 ± 



0.20 5.1 
0.06 15.8 
0.13 7.4 



0.11 
0.11 
0.14 
0.22 
0.19 
0.25 
0.23 
0.26 
0.23 
0.11 



8.9 
9.1 
7.2 
4.5 
5.3 
4.0 
4.4 
4.4 
4.5 
8.9 



Again, as for hydrogen, the many multiplet elements indi- 
cate a significantly distorted dipole pulsation. 



S.S.2 Pulsation properties from abundance patches 

We identified the 'spot' line profile patterns in Fig. [10] as 
A (rotational phases 0.2 — 0.8) and B (rotational phases 
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Figure 11. Amplitude and phase variations for Features A 
(closed circles) and B (open boxes) in Fig. 10 for Nd ill 5102.41 A 
(top) and Nd ill 6145.07 A (bottom), measured with Gaussian pro- 
file fits. A superposed triplet dipole pattern is shown for compar- 
ison with a solid line. The group of measurements near the mag- 
netic crossover phase (^rot = 0.75 occur during the double-lined 
phase and show a phase shift near tt rad between the two poles, 
as expected for a basically dipolar pulsation mode. 



0.7 — 1.3). All spectra were re-measured by fitting Gaussian 
profiles to the lines, while tracing the mean movement of the 
A and B features individually with the rotation phase. For 
time intervals of 1 hr the pulsation amplitude and phase were 
measured for fixed u — 1.5576482 mHz, the actual pulsation 
frequency of the distorted dipole mode. 

The two panels in Fig. [TT] show the resulting amplitude 
and phase variations with rotation phase for the two lines 
Nd III 5102.41 A and Nd ill 6145.07 A. Measurements for each 
spot are indicated with different symbols - open boxes for 
feature B and closed circles for feature A. Superposed on 
the figure is a dipole model, constructed using a triplet fre- 
quency fit to the full Nd series. Several details are notable 
in the figure: the pulsation amplitude never reaches zero for 
any of the spots even when at the stellar limb; the modulated 
amplitude maximum for spot A near rotation phase 0.5 (the 
spectra only cover the ingress part), is sharper, steeper and 
shifted in rotation phase with respect to what is expected 
for a pure dipole model. This is the result of concentration of 



Ndiii towards the magnetic and pulsation poles (see Fig. [8). 
During quadrature at rotation phase 0.75, it is seen that the 
pulsation amplitudes for both spots are similar, and about 
2/3 of the maximum of the amplitude modulation, a conse- 
quence of the projection factor (see below). The increased 
scatter at this phase is due to line blending between the 
double- line components, and weaker absorption lines, as the 
Nd abundance is less overabundant over the visible hemi- 
sphere of the star at this phase. Note however, how clearly 
the pulsation phases divide for the two spots, separated by 
^ TV rad, clearly indicating the global / = 1 dipole pulsation. 

With both spots at the limb of the star, and with the 
ability to study them separately, it is interesting to ask 
whether we can see any evidence of a horizontal component 
to the pulsations. For example, if the pulsation were purely 
horizontal, rather than vertical (radial), we would expect to 
see maximum pulsation amplitude at this rotational phase 
(0.75) for the Ndiii lines. On the other hand, because the 
Ndiii spots are large, there is a component of pulsation in 
the direction of the line-of-sight for purely vertical pulsation. 
A simple integral over the half of the visible hemisphere in 
which a spot is seen (see, e.g., |Parsons 1972 for the simpler 
case of Cepheid pulsation) indicates that for purely vertical 
pulsation velocities we will still expect to see an amplitude 
about 2/3 that seen at pulsation maximum when the spot is 
closest to the line-of-sight. This is what we do see in Fig.fTTl 
hence there is no obvious evidence of a horizontal component 
to the velocity field. We will model this in more detail in a 
future pubUcation to put limits on the horizontal compo- 
nent. However, as mentioned in Sect. 1.3 [Sli3 (l2005l ) shows 
the horizontal motion is expected to be zero at the magnetic 
poles. Thus the best strategy for detecting horizontal motion 
is through chemical spots located away from the poles. 

Figure[T2] shows the line profile variability in 
Nd III 6145.07 A for all 32 1-hr bins that we divided 
the spectra into. Each sub-panel shows the average spec- 
trum at this rotation phase, the residual variability for 
the spectra phased with the pulsation phase progressing 
downwards. The lower insert of each panel is the standard 
deviation for the individual series of spectra, pixel-by-pixel. 
Flux is in all cases in units of normalised flux of the 
stellar continuum. The figures show a tt rad pulsation phase 
reversal from rotation phase iprot = 0.0 to 0.5 (seen as a 
colour inversion) that can be seen to progress with the 
stellar rotation^ 

Recently, IShibahashi et al.l (|2008l l have interpreted 
blue-to-red running wave features in the line profile vari- 
ability in the roAp star 7 Equ as a manifestation of an up- 
per atmospheric shock wave. Similar blue-to-red features are 
also seen in the line profile varia tions of other roAp stars 
(see, e.g., iKochukhov et al.ll2007l ). However, up to now ob- 
servations have generally been made for short time spans of 
typically 2 h and do not cover different rotational phases. An 
important exception to this is the pioneering study of line 
profile variability in the roAp star HR3831 by Kochukhov 
(2006; see his figures 2 and 3) which shows similarities in 
the behaviour to that seen in Fig. [12] for HD 99563. 

We see in Fig. [12] the line profile variations of a spectral 
line simply shifting in wavelength without a change of shape 
at phases near pulsation maximum when the pulsation axis 
is closest to the line-of-sight, as is expected for a dipole pul- 
sation mode at this aspect. This is clear in the top row of 
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Figure 12. Residual line profile variations in Ndlll 6145.07 A for the 32 time series bins, with rotation phase increasing to the right and 
downwards, each subtracted from the average spectrum (shown on top of each panel) of that bin. Pulsation phase increases downwards 
in the middle panels. The standard deviation of each bin series is shown in the bottom of the panels. The colour codes are —0.01 (black) 
to +0.01 (white) in normalised intensity. 



the figure at phase, e.g., (prot = 0.046 and other phases near 
it, and again in the bottom right panel of the figure at phase 
(frot — 0.501. Yet remarkably at magnetic quadrature when 
the Ndiii line is doubled and we can examine the individ- 
ual components seen at an aspect angle near 90°, there is a 
clear blue-to-red running wave component visible. Note par- 
ticularly in Fig. [12] (continued) the upper row of panels and. 



e.g., (/?rot = 0.7 33. This same beh aviour was reported first 
for HR3831 bv iKochukho^ (|2006l ) (who also reported red- 
to-blue features which we possibly see for HD 99563 at some 
rotation phases, e.g. 0.146 in Fig. 12) and is clearly impor- 
tant for our understanding of these running wave features 
seen both here and in other roAp stars. 

Further interpretation of these features for Nd, as well 
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Figure 12. Fig. 12 continued. Note the clear line doubling for the rotation phases in the top row, observed near the time of magnetic 
quadrature. 



as for other lines and for other stars, requires modelling of 
the pulsations, which is work now in progress. Whatever 
the cause of the blue-to-red running wave line profile varia- 
tions - whether shock waves as proposed and modelled by 
IShibahashi et all (|2008l ) , or some other cause - the fact that 
they are aspect dependent with rotation in HD 99563 and 
HR3831 is an important characteristic that needs to be ex- 
plained. 



3.4 Europium 

This element shows a highly inhomogeneous surface distri- 
bution in the Doppler imaging map seen in Fig. 1131 Two 
circumpolar patches are viewed most favourably near the 
line-of-sight around (frot = 0.45 and 0.95; the third spot is 
rather concentrated and appears nearest the line-of-sight at 
(frot = 0.65. 



3.4-1 Observed line profile variations 

In the sequence of spectra in Fig. ll4[ the two strongest Eu 
lines Eu II 6437.64 A and Euii 6645.06 A show three moving 
features which are identified in the figure as A, B and C. 

Tabled shows the results of frequency analyses of three 
elements, Nd, Eu and Ca. The radial velocities were mea- 
sured with Gaussian fitting or the CoG procedure (for Ca) 
to the full line profiles in all spectra. To avoid spectra with 
mixed contributions, some ~40 spectra around the magnetic 
crossover phase (cp = 0.75) were not used. Detrending for low 
frequencies was performed with first and second order poly- 
nomial fitting to spectra dominated by the A or B features 
separately (third and fourth order polynomials for Eu and 
Ca) and up to 3 low frequencies were prewhitened. The re- 
sults indicate that Nd requires a nonuplet solution and thus 
has a high distortion to the oblique dipole modulation, while 
Eu has little power left after a triplet fit. Ca may, however, 
require at least a quintuplet pattern , but as it has a very low 
amplitude of a few hundred km s~^ (|Elkin, Kurtz Mathvsl 
|2005|), the frequencies are not significant when fitting the 
full nonuplet solution. These results are consistent with the 
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Figure 13. Doppler imaging map for Eu ll 6645.06 A. The abun- 
dances are on a scale of log Nh = 12; the solar abundance of 
Eu is 0.48 on this scale (Asplund et al. 2005). Thus the highest 
overabundance of Eull in the spots rises to 10^ times solar at 
the centres of the three spots, two of which are offset from the 
magnetic and pulsation poles, and the third which is isolated at 
longitude 290° . The white colour below latitude —44° represents 
the part of the star that is never visible because of the i = 44° 
rotational inclination. 

abundance distributions of these elements. The small A and 
B spots of Eu do not sample much of the visible hemisphere 
of the star, except near the poles, so they are fitted well by a 
simple dipole model. That means that the distortions to the 
dipolar pulsation that give rise to the other members of the 
nonuplet for Nd come from regions further from the poles 
than the Eu samples. These differences may eventually be 
exploited for more detailed models of the pulsation velocity 
field. The high 5kms~^ amplitude seen for Eu is consistent 
with its concentration towards the poles, and may also in- 
dicate that it is radiatively stratified to high levels in the 
atmosphere where the pulsation amplitude is higher. The 
phase differences (relative to Eq.|4| of the two main multi- 
plet sidelobes of Eu in Tableware, such as shown for Nd, at 
least partly due to different surface distributions than H. 

3.4-2 Pulsation properties from abundance patches 

Angular phase- amplitude variation was studied for the Eu 
A and B features in the same manner as done for Nd. The 
phase differences between Nd and Eu of the A and B fea- 
tures in Figs. [11] and [15] indicate different surface distribu- 
tions such that different parts of the pulsation geometry are 
mapped, consistent with the Doppler imaging maps. 

The Eu C spot was also studied for angular phase- 
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Figure 14. Line region in the vicinities of Eull 6437.64 A (left) 
and Eull 6645.06 A (right) shown for all spectra sorted according 
to rotation phase of HD 99563. Absorption features appear darker 
and in black. Four phase intervals were not covered by fast spec- 
troscopy and are shown as smooth dark gray bands. Two line 
features relate to dominant surface patches of high Eu concen- 
tration. These are referred to as A (rotational phases 0.2 — 0.8) 
and B (0.7 — 1.3). A third, much more shallow feature is seen at 
rotational phases 0.7 — 0.8. 

amplitude variation, but the amplitude was found low and 
rather constant. This is not surprising, since the C spot (see 
Fig.[l3| lies about 40° away from the longitude of the pul- 
sation pole where significantly lower pulsation amplitude is 
expected for a dipole mode. 



4 CONCLUSIONS 

We have performed a first analysis of the pulsation of 
HD 99563 based on one of the most comprehensive spectro- 
scopic efforts on a roAp star so far. We show that the pul- 
sation is well-described by a single distorted oblique dipole 
mode, with different distortion terms for different elements. 
In the future, we will use these differences to make further 
inference concerning the pulsation velocity field - ultimately 
in 3-D. By using Doppler imaging maps of the surface distri- 
bution of Nd and Eu, we have identified two relatively con- 
centrated spots near the magnetic poles. During rotation 
phases at magnetic crossover, line doubling in absorption 
features of these elements (with each component relating to 
a spot) is considerable, and we show that the stellar photo- 
sphere at opposite sides of the star have mutual tt rad offset, 
again supporting the oblique dipole rotator geometry. These 
spots allow the study of the pulsation mode for only parts of 
the visible hemisphere of the star, and from varying aspect, 
an opportunity only possible for roAp stars. Remarkably, 
blue-to-red running wave features were found only near the 
rotation phases where the lines were doubled and the spots 
were seen at high aspect angle. 

Our study included a short data set we also obtained 
with UVES on the VLT 3 years earlier (Elkin et al. 2005), 
which after a minor revision of the pulsation frequency 
showed amplitudes and phases for H, Nd and Eu that in 
general are consistent with the new data. Some roAp stars 
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Table 4. Pulsation amplitudes and phases (ip) fitted to frequencies calculated relatively to a slightly revised u and 
to Urot in TableO The data were detrended for low frequency trends and low frequencies. The radial velocities were 
measured using Gaussian fitting, except for lines of calcium where the CoG procedure was used. Note the high 
amplitude of the Eu lines. 



Id. 


Amplitude 




SjN Amplitude 




SjN Amplitude 




SjN 




(ms-l) 


(rad) 




(ms-l) 


(rad) 




(ms-l) 


(rad) 






Ca 1 6439.07 A 




Eul 


16645.06 A 




Eull 


6437.64 A 




u 


190 ± 139 


-1.39 ±0.72 


1.4 


452 ±419 


-1.65 ±0.91 


1.1 


264 ± 585 


1.46 ±2.21 


0.5 




409 ± 143 


1.46 ±0.35 


2.9 


2631 ± 466 


-1.59 ±0.17 


5.6 


2292 ± 656 


-1.18 ±0.28 


3.5 




342 ± 144 


2.90 ±0.42 


2.4 


1996 ± 453 


-1.09 ±0.23 


4.4 


1931 ± 646 


-1.41 ±0.34 


3.0 


I/ — 2z/rot 


240 ± 100 


1.97 ±0.41 


2.4 


377 ± 320 


0.29 ±0.85 


1.2 


488 ± 447 


-2.58 ±0.92 


1.1 




120 ± 99 


1.36 ±0.82 


1.2 


292 ± 321 


0.66 ± 1.10 


0.9 


941 ± 447 


0.10 ±0.48 


2.1 


ly — SUrot 


263 ± 91 


2.62 ±0.34 


2.9 


289 ± 275 


-0.10 ±0.94 


1.0 


271 ± 385 


0.69 ± 1.42 


0.7 




95 ±91 


-0.58 ±0.95 


1.1 


215 ± 274 


2.90 ± 1.27 


0.8 


784 ± 384 


2.38 ±0.49 


2.0 




203 ± 88 


-2.24 ±0.44 


2.3 


802 ± 278 


2.01 ±0.34 


2.9 


120 ± 383 


-1.99 ±3.22 


0.3 




11 ±88 


3.14 ±8.28 


0.1 


422 ± 276 


1.25 ±0.66 


1.5 


689 ± 386 


-2.62 ±0.55 


1.8 


2v 


95 ±37 


-0.10 ±0.39 


2.6 


141 ± 114 


-2.85 ±0.81 


1.2 


279 ± 160 


0.33 ±0.58 


1.7 


2i/-2l/rot 


40 ±37 


-0.10 ±0.92 


1.1 


299 ± 114 


1.25 ±0.38 


2.6 


46 ± 161 


2.34 ±3.50 


0.3 




Nd III 5102.41 A 




Nd III 6145.07 A 




Ndiii 6327.24 A 




V 


281 ± 68 


-0.29 ±0.24 


4.2 


367 ± 186 


-0.95 ±0.51 


2.0 


267 ± 83 


-0.99 ±0.31 


3.2 


i/-l/rot 


1000 ± 78 


2.78 ±0.08 


12.7 


1339 ± 193 


-2.89 ±0.14 


6.9 


1439 ± 86 


-2.73 ±0.06 


16.7 




530 ± 78 


2.64 ±0.15 


6.8 


778 ± 196 


1.72 ±0.25 


4.0 


282 ± 89 


1.52 ±0.32 


3.2 


l/-2l/rot 


318 ±57 


1.16 ±0.18 


5.5 


597 ± 132 


0.33 ±0.22 


4.5 


456 ± 60 


0.42 ±0.13 


7.6 




472 ± 57 


-2.43 ±0.12 


8.2 


974 ± 134 


-1.70 ±0.14 


7.3 


480 ± 59 


-1.37 ±0.12 


8.1 


Z/-3z/rot 


194 ± 47 


-1.40 ±0.24 


4.2 


339 ± 121 


-0.89 ±0.36 


2.8 


381 ± 54 


-1.44 ±0.14 


7.1 


l/+3z/rot 


242 ± 47 


-1.40 ±0.19 


5.2 


450 ± 123 


-0.67 ±0.27 


3.7 


286 ± 54 


-0.91 ±0.19 


5.3 


V — Av^ot 


25 ±50 


0.78 ±1.97 


0.5 


305 ± 117 


2.49 ±0.38 


2.6 


346 ± 52 


2.84 ±0.15 


6.6 


l/+4z/rot 


208 ± 50 


0.66 ±0.24 


4.2 


575 ± 121 


1.75 ±0.21 


4.8 


340 ± 55 


1.96 ±0.16 


6.2 


2v 


146 ± 19 


1.05 ±0.13 


7.6 


113 ±49 


0.33 ±0.44 


2.3 


59 ±22 


-0.04 ±0.37 


2.7 




147 ± 19 


0.83 ±0.13 


7.6 


144 ± 49 


2.77 ±0.34 


2.9 


160 ± 22 


1.39 ±0.14 


7.3 



exhibit short lifetimes of pulsation modes - in the most ex- 
treme case, HD 60435, the mode lifetimes are less than a 
week. The stability of HD 99563 over a 3-yr time span sim- 
plifies the interpretation of the pulsations. This also suggests 
that it will be possible to merge future observations of the 
rotational phases not covered by our data for a more detailed 
study. 

We found a clear rotational variation in the radial veloc- 
ities of the Ha line in addition to the pulsational variations 
which we interpret to be the first observational evidence of 
hydrogen overabundance at the magnetic/pulsational poles, 
as expected theoretically as a result of helium depletion. 
While we prefer this interpretation, alternatively, we found 
the first planet orbiting a roAp star, a hot Jupiter with a 
mass of 26 times the mass of Jupiter. In future work we will 
model the Ha line radial velocity variations in detail to test 
our preferred hypothesis of hydrogen spots. 
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Figure 15. Amplitude and phase variations for features A 
(closed circles) and B (open squares) for Eull 6437.64 A (top) and 
Eull 6645.06 A (bottom). The line profile fitting was performed 
using the CoG procedure which gives a lower amplitude than with 
Gaussian fitting (such as in Fig. ll2|) . but the amplitude exceeds 
that of Ha. The superposed triplet dipole model was fitted to the 
radial velocities of the strongest Eu line, Eull 6645.06 A. The fea- 
ture C (not shown) was also studied, but had very low amplitude 
and a high scatter in phases (consistent with its longitude) and 
has been excluded in this figure for clarity. 
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